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ABSTRACT 

We present moderately-high resolution echellc observations of the nucleus of 
NGC 5195, the line of sight to which samples intervening interstellar material 
associated with the outer spiral arm of M51. Our spectra reveal absorption from 
interstellar Na I, K I, Ca n, and CH + , and from a number of diffuse interstellar 
bands (DIBs), at a velocity close to that exhibited by H I 21 cm emission from 
M51 at the position of NGC 5195. The H I column density implied by the equiv¬ 
alent width of the A5780.5 DIB, based on the relationship between 11^(5780.5) 
and N (H i) derived for sight lines in the local Galactic interstellar medium, is 
consistent with the column density obtained from the integrated H I emission. 

The H 2 column density predicted from the observed column density of K I, using 
the Galactic relationship between N (K 1 ) and iV(H 2 ), is comparable to N (H 1 ), 
suggesting a high molecular fraction (~0.65) for the M51 gas toward NGC 5195. 

The DIBs toward NGC 5195 are, on average, ~40% weaker than would be ex¬ 
pected based on the K I column density, a further indication that the gas in 
this direction has a high molecular content. The M51 material is characterized 
also by a high A(Na l)/A(Ca 11 ) ratio (> 11), indicative of a high degree of Ca 
depletion, and a high R / a(5797.1)/II / a( 5780.5) ratio (~1.6), suggestive of either 
a very weak ambient radiation held or a significantly shielded environment. A 
high A(CH + )/A(CH) ratio (> 2.3) for the M51 material toward NGC 5195 may 
be the result of enhanced turbulence due to interactions between M51 and its 
companion. 
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1. INTRODUCTION 

The grand design spiral galaxy M51 (NGC 5194) has long been known to be interacting 
with its companion, the SBO galaxy NGC 5195, which is seen in projection behind dusty 
interstellar gas now understood to be associated with the northeastern spiral arm of M51. 
Spinrad (1973) commented on the obvious dust lanes that are seen in front of NGC 5195, 
and speculated that some of these could be due to foreground absorption from the outer 
arm of M51. However, it was Van Dyk (1987) who confirmed this speculation, discovering 
anomalously strong Na I D lines attributable to interstellar gas in the arm of M51 projected 
onto the nucleus of NGC 5195. The strength of the interstellar features combined with 
the brightness of the background source makes the line of sight to NGC 5195 ideal for 
studying the interstellar medium (ISM) of a nearby (d ~ 8.0 Mpc; Karachentsev et al. 2004) 
interacting galaxy using the technique of absorption-line spectroscopy. The low resolution 
of the spectrum analyzed by Van Dyk (1987), however, precluded a detailed investigation of 
the ISM in this direction. 

In Ritchey & Wallerstein (2012), we presented moderately-high resolution echelle obser¬ 
vations of SN 201 ldh in M51, which revealed numerous interstellar absorption components 
in Na I D and Ca II H and K associated with the supernova host galaxy. The M51 com¬ 
ponents were found to span a velocity range of over 140 km s' 1 , making them similar to 
intermediate-velocity clouds observed in the halo of the Milky Way. Yet many of the M51 
components toward SN 201 ldh, even those at extreme positive and negative velocities, were 
found to have Na i/Ca II column density ratios of ~1.0. Such values are typically associated 
with relatively cool quiescent gas in the disk of a galaxy like M51 or the Milky Way. Some 
of these components may thus be tracing gas that has been tidally stripped from the disk 
of M51 or its companion (or both) during one of the multiple encounters which are thought 
to have occurred recently between the two galaxies (e.g., Salo & Laurikainen 2000; Theis & 
Spinneker 2003). 

Background sources suitable for probing the ISM of nearby galaxies that are bright 
enough for high-resolution absorption-line spectroscopy are rare. Extragalactic studies using 
individual stars as background sources are currently limited to examining sight lines within 
Local Group galaxies, such as the Large and Small Magellanic Clouds (LMC and SMC; e.g., 
Cox et al. 2006, 2007; Welty et al. 2006; Welty & Crowther 2010), M31 (e.g., Cordiner et 
al. 2008a, 2011), and M33 (e.g., Cordiner et al. 2008b; Welsh et al. 2009). Bright supernovae 
are perhaps the best sources for probing the ISM of galaxies beyond the Local Group, 
and, indeed, many nearby galaxies have been investigated through supernova absorption¬ 
line spectroscopy (e.g., Wallerstein et al. 1972; D’Odorico et al. 1989; Bowen et al. 1994, 
2000; Vladilo et al. 1994; Ho & Filippenko 1995; Sollerman et al. 2005; Cox & Patat 2008, 
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2014; Welty et al. 2014; Ritchey et al. 2015). However, supernovae fade over time, making 
follow-up high-resolution observations difficult, and the occurrence of any single supernova 
cannot be predicted. The nuclei of nearby spiral galaxies can be relatively bright, but often 
do not show measureable interstellar absorption features (unless the galaxy is observed in an 
edge-on orientation). Thus, the chance position of NGC 5195 located behind a prominant 
spiral arm of M51 presents a unique situation. 

In this paper, we examine moderately-high resolution echcllc observations of the nucleus 
of NGC 5195, which allow us to study the properties of the interstellar gas seen in absorption 
along the line of sight. The spectra reveal redshifted interstellar absorption features of Na I, 
K I, Ca II, and CH + , along with a number of diffuse interstellar bands (DIBs), all of which 
help to constrain the physical conditions in the intervening interstellar material in front of 
NGC 5195. We briefly examine some of the prominent emission lines and stellar absorption 
lines in the spectrum of NGC 5195, and in spectra of the nucleus of M51, in order to 
place our measurements of the interstellar features in the proper context. In particular, a 
comparison of the radial velocities of the various absorption and emission features allows 
us to unambiguously associate the redshifted interstellar absorption lines toward NGC 5195 
with gas in the spiral arm of M51. 


2. OBSERVATIONS AND DATA REDUCTION 

Our observations were acquired using the Astrophysical Research Consortium (ARC) 
echelle spectrograph on the 3.5 m telescope at Apache Point Observatory (APO). The echcllc 
provides complete wavelength coverage in the range 3800-10200 A with a resolving power of 
R 31,500 (Av ~ 9.5 km s -1 ). A total of eighteen 30 minute exposures of the nucleus of 
NGC 5195 (V = 9.55 ± 0.07; Gil de Paz et al. 2007) were obtained between 2012 January 
and 2013 April. In the earliest spectra that we obtained, it was apparent that the Na I D 
features at the redshift of M51 were quite strong in this direction, while the Ca II H and K 
lines were buried in the noise. Since the Na i/Ca II column density ratio can be a useful 
diagnostic of physical conditions, we continued to accumulate observations to build up the 
signal until the Ca II lines were clearly visible. On most occasions, the nearby, early-type 
star r) UMa (B3 V; V = 1.85) was observed on the same night as NGC 5195 to aid in the 
removal of telluric absorption lines, which affect the Na I D and K I AA7664, 7698 features. 
We obtained two 30 minute exposures of the nucleus of M51 (V = 8.36 ± 0.06; Gil de Paz et 
al. 2007) in 2012 January for comparison with NGC 5195. 

Reduction of the raw data was accomplished using a semi-automated pipeline reduc- 
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tion script, which employs standard IRAFhl routines for bias correction, cosmic-ray removal, 
scattered-light subtraction, one-dimensional spectral extraction, flat-fielding, and wavelength 
calibration. The atmospheric absorption lines near the Na I D and K I features were removed 
from the individual, calibrated exposures using the TELLURIC task within IRAF. We find 
that this procedure is effective even for removing the strong absorption lines associated with 
the A band of atmospheric O 2 in the vicinity of the K I doublet, which normally hinder the 
analysis of the K 1 A7664 line. For these observations, we are aided also by the redshift of 
M51, which places the associated K 1 A7664 feature in a region mostly free from absorption 
by atmospheric O 2 . After correcting for telluric absorption lines, all of the exposures were 
shifted to the reference frame of the local standard of rest (LSR) and co-added to produce 
a final, high signal-to-noise ratio (S/N) spectrum. Any overlapping portions of the echcllc 
orders containing absorption (or emission) lines of interest were combined to improve the 
S/N and the resulting spectra were normalized to the continuum via low-order polynomial 
fits. 


The S/N (per pixel) in the final, co-added spectrum of NGC 5195 peaks at ~200 near 
the K I doublet, and decreases to ~150 near the Na I D features. At the Ca II H and K lines, 
the S/N is only ~15, while near the Ca II infrared triplet, the S/N is ~115. For comparison, 
the S/N in our spectrum of the nucleus of M51, which had a total integration time of only 60 
minutes, peaks at ~60 near K I and decreases to ~10 at Ca II K. Figure 1 presents a portion 
of the final, co-added spectrum of NGC 5195 in the vicinity of Na I D. In this figure, the very 
narrow absorption features near the rest wavelengths of the D lines are due to interstellar 
Na I in the Milky Way. The broader, redshifted absorption features, which reach almost zero 
intensity at the core, are presumably due to interstellar material in the spiral arm of M51, 
which is seen in projection against NGC 5195. The M51 components are superimposed onto 
very broad stellar Na I D lines, which are further redshifted and presumably arise from red 
giant stars in the nucleus of NGC 5195. 

In addition to Na I D, we clearly detect the interstellar K I and Ca II doublets to¬ 
ward NGC 5195 near the velocity expected for gas associated with M51 (i.e., near ulsr ~ 
430 km s -1 ; see Section 3.1). The narrow Galactic component seen in Na I D absorption 
(Figure 1), which has a velocity of n LSR ~ —30 km s -1 , is also detected in K I and absorp¬ 
tion near this velocity can be seen in Ca II K, as well, ffowever, the main focus of this 
investigation will be on the interstellar absorption features that we associate with the spiral 
arm of M51. We detect a strong line of interstellar CH + A4232 toward NGC 5195 near the 


2 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso¬ 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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expected M51 velocity, but do not see evidence for CH A4300. (The non-detection of CH 
is not that unexpected given that the S/N in this region is only ~30.) We do not detect 
any interstellar absorption features toward the nucleus of M51 other than a narrow Na I D 
component associated with the Milky Way, which is weaker than the corresponding feature 
toward NGC 5195. 

Eight DIBs are detected toward NGC 5195 at the same velocity as the other interstellar 
species that we presume to be associated with M51. The detected DIBs are those with rest 
wavelengths of 5780.5, 5797.1, 5849.8, 6196.0, 6203.6, 6376.1, 6379.3, and 6613.6 A (where 
we have adopted the convention of Friedman et al. (2011) of quoting DIB wavelengths to one 
decimal place). Searches for other prominent DIBs in the NGC 5195 spectrum, such as the 
others studied by Friedman et al. (2011) or those detected by Sollerman et al. (2005) toward 
SN 2001el in NGC 1448, were unsuccessful. However, relatively weak and/or significantly 
broadened DIBs would be difficult to positively identify in our moderate S/N data, given the 
shape of the underlying continuum, which comprises many broad and often blended stellar 
absorption features. The DIBs that we do identify are generally much narrower than the 
underlying stellar features, as illustrated in Figure 2, which shows a portion of the NGC 5195 
spectrum in the vicinity of the (redshifted) A5780.5 and A5797.1 DIBs. Also plotted in Fig¬ 
ure 2 is a smoothed spectrum of the prototypical red giant star Arcturus (obtained from the 
atlas spectrum of Arcturus published in electronic form by Hinkle et al. 2000), which provides 
a template for the underlying stellar absorption features in the NGC 5195 spectrum. The 
close correspondence between the smoothed Arcturus spectrum and our adopted continuum 
fit for the region surrounding the A5780.5 and A5797.1 DIBs indicates that the continuum 
has been reasonably well fit and that any contamination of the DIB absorption by stellar 
lines is minimal. Similar comparisons for the regions surrounding the other detected DIBs 
indicate that the continua in those regions have also been fit appropriately. 

In Figure 3, we present continuum-normalized spectra for the interstellar atomic and 
molecular species observed in the direction of NGC 5195. We also show (in the top panel 
of Figure 3) the H I 21 cm emission profile of M51 at the location of NGC 5195 provided 
by The H I Nearby Galaxy Survey (THINGS; Walter et al. 2008). This spectrum was 
extracted from the THINGS data cubto for M51 using the average of the four pixels nearest 
to the coordinates of NGC 5195. Figure 4 presents the absorption profiles for the eight 
DIBs identified in our co-added spectrum of NGC 5195. In Figure 5, we show continuum- 
normalized emission profiles for the [N u] AA6548, 6583, Ha, and [S ii] AA6716, 6730 emission 
lines, which are prominent in the NGC 5195 spectrum. In Figure 6, we present the absorption 


3 Available from the THINGS website at: http://www.mpia-hd.mpg.de/THINGS/ 




profiles of the Ca II infrared triplet, which (presumably) traces the red giant population in 
NGC 5195. 


3. ANALYSIS 

3.1. Radial Velocities and Line Widths 

Radial velocities and line widths for the prominent emission lines and for several of 
the stellar absorption lines seen in the spectrum of NGC 5195 and of the nucleus of M51 
were determined by fitting Gaussian functions to the observed features within IRAF (see 
Tables 1 and 2). For both galaxies, we note a relatively small, yet statistically significant, 
offset between the radial velocities of the emission lines and those of the stellar absorption 
features. For NGC 5195, the mean centroid velocity of the stellar absorption lines listed in 
Table 1 is 616 ± 9 km s' 1 , while the mean velocity of the emission lines is 570 ± 7 km s' 1 . 
(Errors for mean values reported in this paper are la standard deviations.) The overall 
mean for the emission lines and stellar absorption lines together is 588 ± 25 km s' 1 , which 
is in good agreement with the velocity listed for NGC 5195 in the Updated Zwicky Catalog 
(UZC): ulsr = 584 ±43 km s' 1 (Falco et al. 1999, applying a heliocentric-to-LSR correction 
of 12 km s' 1 ). In the case of M51, the emission lines listed in Table 2 show a mean velocity 
of 502 ± 7 km s' 1 , while the stellar absorption lines have a mean of 477 ± 7 km s' 1 . The 
overall mean is 495 ± 14 km s' 1 , which (likewise) is consistent with the UZC value for M51: 
ulsr = 477 ± 61 km s' 1 . 

The central velocity of the global H I 21 cm emission profile of M51 from THINGS is 
ulsr = 468 km s' 1 (Walter et al. 2008, applying the same heliocentric-to-LSR correction 
as above). At the position of NGC 5195, however, the velocity held (moment 1 map) of 
M51 provided by THINGS indicates an intensity-weighted mean H I velocity of ulsr — 
434 km s' 1 . The mean column-density weighted velocity of the interstellar atomic and 
molecular absorption lines seen in the direction of NGC 5195 is 411±4 km s' 1 , and (similarly) 
the mean velocity exhibited by the DIBs is 411 ± 5 km s' 1 (see Tables 3 and 4). Since these 
velocities are close to the mean H I velocity from THINGS at the position of NGC 5195, 
we can unambiguously associate the interstellar features with gas in M51 rather than in 
its companion. We do note, however, the slight difference between the H I velocity and the 
velocities of the interstellar absorption features toward NGC 5195, which is clearly discernible 
in Figure 3. We discuss this apparent discrepancy further in Section 4. 

From the widths of the stellar absorption lines and the emission lines listed in Tables 
1 and 2, we can derive estimates for the stellar velocity dispersion (cr*) and the gaseous 
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velocity dispersion ( a g ) in NGC 5195 and in the nucleus of M51. For a Gaussian line profile, 
the full width at half maximum (FWHM) is related to the velocity dispersion according 
to FWHM = 2^2 In 2a a; 2.35a. The values of a derived using this relation are listed 
in the last columns of Tables 1 and 2. In converting the observed FWHM to a, we have 
taken into account both the instrumental broadening (which, for our observations, is almost 
negligible), and the intrinsic widths of the stellar absorption lines. Estimates for the intrinsic 
widths of the stellar lines were obtained from the (unsmoothed) atlas spectrum of Arcturus 
(Hinkle et al. 2000), which we assume to be typical of the red giants in the nuclei of these 
galaxies. For NGC 5195, the stellar absorption lines yield a mean velocity dispersion of 
a* = 112 ± 13 km s _1 , while the emission lines show a mean of a g = 70 ± 10 km s -1 . For the 
nucleus of M51, the mean stellar velocity dispersion is a* = 91 ± 16 km s -1 , while the mean 
gaseous velocity dispersion is a g = 72 ± 13 km s -1 . Our estimates of a* for these galaxies 
are consistent with the values adopted by Ho et al. (2009), who give a* = 124.8 ±8.1 km s^ 1 
for NGC 5195 and a* = 96.0 ± 8.7 km s _1 for M51. However, our results also suggest 
that a* > <j g for both NGC 5195 and M51 (although the difference for M51 may not be 
significant). When considered together, the differences between the radial velocities and 
velocity dispersions inferred from the stellar absorption lines and those derived from the 
emission features may be indicative of real differences in the kinematics of the stars and the 
ionized gas in the central regions of these galaxies. 


3.2. Equivalent Widths and Column Densities of the Interstellar Features 

The equivalent widths (W\) of the interstellar atomic and molecular absorption lines 
observed in the direction of NGC 5195 (that are associated with interstellar gas in M51) 
were computed by direct integration over the absorption profiles (see Table 3). The limits of 
integration were set according to the locations where the optical depth reached zero. For dou¬ 
blets, the limits were determined based on the stronger line and applied in an identical way to 
the weaker line. Uncertainties in the integrated equivalent widths include contributions from 
photon noise and continuum placement. The photon noise uncertainties were derived using 
the formula given in Jenkins et al. (1973), which can be written a^ A = (AA)(iV p ) 1 / 2 (RMS), 
where AA is the difference in wavelength between adjacent pixels, N p is the number of pixels 
sampling the absorption profile, and RMS is the root mean square of the noise in the con¬ 
tinuum. The DIB equivalent widths and uncertainties were computed in the same manner 
(see Table 4). As an aid in determining the limits of integration for the DIB equivalent 
widths, we compared the profiles of the DIBs detected toward NGC 5195 with those seen 
toward Galactic stars (e.g., Friedman et al. 2011). The adopted limits are explicitly shown 
in Figure 4. 



Column densities for the atomic and molecular species were determined by means of 
the apparent optical depth (AOD) method, which is a useful approach when it is suspected 
that the absorption profiles may be affected by unresolved saturation (see Savage & Sembach 
1991). While there is evidence for component structure in the atomic and molecular absorp¬ 
tion profiles, in the form of a weak wing redward of the main absorption components (see 
Figure 3), the moderate resolution of our observations and the relatively low S/N achieved 
below 5000 A suggests that a detailed profile synthesis analysis would not be an ideal ap¬ 
proach for these data. However, we can still recover measurements of total column densities 
through use of the AOD method. Apparent column densities (N a ) were derived by integrat¬ 
ing the AOD profiles over the same intervals used to compute the equivalent widths. Values 
of N a for the individual lines detected toward NGC 5195 are listed in Table 3. Uncertainties 
in these quantities were determined through standard error propagation (see, e.g., Prochaska 
et al. 2001), and also include continuum placement uncertainties. 

A comparison between the apparent column densities of the two K I transitions indicates 
that these features may be very mildly saturated. The difference in N a for the two lines 
suggests that a correction of +0.04 dex (applied to the weaker line) is appropriate, according 
to the correction scheme devised by Savage & Sembach (1991). This correction has been 
used to derive the adopted column density (V a d op t) for K I listed in Table 3. A similar 
comparison for the two Ca II transitions indicates that a correction of +0.12 dex is needed, 
and this has been applied to the weaker line to derive the adopted Ca II column density. 
The interstellar Na I lines toward NGC 5195 (that are associated with M51) are very heavily 
saturated. The difference in N a for the two transitions suggests that a correction of at least 
+0.49 dex is necessary. Yet, while this correction has been applied to the weaker line, we 
still regard the resulting Na I column density as a lower limit. No saturation correction is 
needed in the case of CH + A4232. A tentative detection of the weaker CH + A3957 line yields 
a CH + column density consistent with that given in Table 3. 


4. DISCUSSION 

The column densities of the interstellar atomic and molecular species and the equivalent 
widths of the DIBs detected toward NGC 5195 provide a means to estimate some of the basic 
properties of the intervening interstellar material in the spiral arm of M51. Several well- 
defined correlations are known to exist between various interstellar quantities for sight lines in 
the local Galactic ISM (e.g., Welty et al. 2014), and, here, we explore what these correlations 
would imply for the M51 material toward NGC 5195. Specifically, we consider the known 
correlations between the equivalent width of the A5780.5 DIB and the column density of 
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atomic hydrogen iV(H i) (e.g., Herbig 1993; Welty et al. 2006; Friedman et al. 2011), between 
the K I column density and the total hydrogen column density iV(H tot ) = N( Hi) + 21V(H 2 ) 
and/or the H 2 column density (e.g., Welty & Hobbs 2001), between the column densities of 
CH and H 2 (e.g., Welty et al. 2006; Sheffer et al. 2008), and between the equivalent widths 
of various DIBs and E(B — V) (e.g., Friedman et al. 2011; Vos et al. 2011). 

To derive specific estimates for N( H i), 7V(H 2 ), fV(H tot ), and E(B — V ) from the quanti¬ 
ties we measure toward NGC 5195, we use the recently updated Galactic relationships involv¬ 
ing these quantities discussed in Welty et al. (2014) and Ritchey et al. (2015). Our measure¬ 
ment ofW A (5780.5) = 71 ±10 mA implies an H I column density of log N (H i) = 20.75±0.12. 
This value (perhaps fortuitously) is identical to that obtained from the THINGS integrated 
H I (moment 0) map of M51 at the position of NGC 5195; the average of the four pixels 
nearest to the coordinates of NGC 5195 in the THINGS map indicates a column density 
of log N (H i) = 20.75 ± 0.06. In contrast, the total hydrogen column density implied by 
the column density of K I toward NGC 5195 is significantly higher; our measurement of 
log iV(K i) = 12.11 ± 0.01 implies log V(H to t) = 21.43 ± 0.15. Taken together, these esti¬ 
mates for N( H i) and lV(H tot ) would imply an H 2 column density of log 1V(H 2 ) = 21.03±0.18, 
and a molecular fraction of /(H 2 ) = 2fV(H 2 )/fV(H tot ) = 0.79 ± 0.18. A somewhat lower H 2 
column density is obtained if we estimate N( H 2 ) directly from N(K i); in this case, we find 
log N( H 2 ) = 20.73 ± 0.17, which, when combined with our estimate for N( H i), would give 
log fV(H tot ) = 21.21 ± 0.12 and /(H 2 ) = 0.65 ± 0.19. We note that a relatively high value of 
V(H 2 ) is not ruled out by the non-detection of CH, as our 3a upper limit of log 1V(CH) < 13.6 
implies log 1V(H 2 ) < 21.0. However, this upper limit does tend to favor the lower H 2 column 
density estimated directly from N (K i). 

There is reason to suspect that the total hydrogen column density predicted from N( K i) 
toward NGC 5195 may be an overestimate. Deviations from the usual Galactic relationship 
between N( K i) and N( H tot ) have been noted for particular regions in the local Galactic 
ISM and for certain extragalactic sight lines. Many sight lines in the Sco-Oph and Orion 
Trapezium regions, for example, exhibit lower values of N( K i) than would be expected based 
on their total hydrogen column densities, most likely as a result of enhanced radiation fields 
in those environments (Welty & Hobbs 2001). Low values of N(K i), relative to V(H tot ), are 
also observed for sight lines in the LMC and SMC; in these cases, the lower values probably 
reflect the combined effects of stronger radiation fields and lower overall metallicities (Welty 
& Hobbs 2001; Welty 2014). Differences in metallicity and/or in the strength of the radiation 
held for the M51 material toward NGC 5195 could therefore (potentially) be affecting our 
estimate of V(H tot ) from N(K i). The H 2 column density predicted from N( K i) should not 
be as affected, however, since there seems to be less regional variation in the N( K i) versus 
N( H 2 ) relation, probably because both K I and H 2 can be affected by the local radiation 
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field (D. Welty 2014, private communication). We note also that a solar metallicity for the 
M51 material is suggested by the fact that the equivalent width of the A5780.5 DIB, when 
compared with the H I column density derived from THINGS data, is consistent with the 
Galactic relationship between Ha( 5780.5) and N (H i) (e.g., Welty et al. 2006; Friedman et 
al. 2011). In the discussion that follows, we will adopt the value of iV(H 2 ) estimated directly 
from IV(K i) (i.e., log IV(H 2 ) ~ 20.73), along with the corresponding values for IV(H to t) and 
/(H 2 ) assuming log IV(H i) ~ 20.75. 

An estimate for the interstellar reddening associated with the M51 material toward 
NGC 5195 can be obtained from the strengths of the various DIBs observed in that direction. 
Again, we use the updated correlations between Ba(DIB) and E(B — V ) reported in Welty 
et al. (2014), who studied all of the DIBs we detect toward NGC 5195 except for A5849.8 
and A6376.1. The measured equivalent widths of the other six DIBs yield values of E(B — V ) 
in the range 0.13-0.60 mag, with a mean of 0.35 ± 0.20 mag (see Table 4). However, in 
Galactic sight lines with high molecular fractions, many of the typically stronger, well- 
studied “standard” DIBs (such as those detected toward NGC 5195) tend to be weaker than 
expected compared to the general trends with IV(H tot ) and E(B — V), but are not weaker 
compared to trends with IV(H i) (Welty et al. 2014), suggesting that the standard DIBs 
trace primarily atomic gas. This may help to explain the weakness of the A5780.5, A6196.0, 
A6379.3, and A6613.6 DIBs toward NGC 5195, which are weaker by 75%, 70%, 52%, and 
44%, respectively, compared to the strengths predicted from the amount of K I present (see 
Table 4). The strengths of the A5797.1 and A6203.6 DIBs toward NGC 5195 are consistent 
with the predictions based on IV(K i), considering the uncertainties. Still, the DIBs toward 
NGC 5195 are, on average, about 40% weaker than predicted, meaning that the E(B — V ) 
values derived from the DIB equivalent widths are likely to be underestimates. 

We can refine our estimate for E(B — V) by considering the relationships that exist 
between the residuals of the observed values of HA (DIB) with respect to the mean Galactic 
trends with both IV(K i) and E(B — V ) (Welty et al. 2014)0 All of the DIBs studied by Welty 
et al. (2014) exhibit positive slopes in (logarithmic) plots of residual with respect to E(B — V) 
versus residual with respect to IV(K i), indicating that when a given DIB is weak relative to 
IV(K i), it will also be weak relative to E(B — V). Use of the relationships for the residuals 
provided by Welty et al. (2014) for the DIBs detected toward NGC 5195 yields the values 
of E(B — V) given in the last column of Table 4. Compared to the original predicted values, 
these adjusted values exhibit a narrower range (0.35-0.55 mag), along with a higher mean 


4 For a given DIB, the “residual” is here defined as the difference between (the logarithm of) the mea¬ 
sured equivalent width and (the logarithm of) the equivalent width predicted by the appropriate Galactic 
relationship involving either IV(K i) or E(B—V). 
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value and smaller standard deviation (0.46 ± 0.09 mag), and should be more representative 
of the reddening toward NGC 5195 that is associated with intervening interstellar material 
in M51. 

Our estimates of log A^(H to t) ~ 21.21 and E(B — V) ~ 0.46 mag for the M51 material 
toward NGC 5195 yield a gas-to-dust ratio of N(H tot )/E(B — V) ~ 3.6 x 10 21 cm -2 mag -1 , 
somewhat lower than the standard Galactic gas-to-dust ratio of 5.6-5.9 x 10 21 cm -2 mag -1 
(Bohlin et al. 1978; Rachford et al. 2009; Welty et al. 2012). For the atomic gas alone, we find 
a gas-to-dust ratio of 1V(H I )/E(B — V) wl.2x 10 21 cm -2 mag -1 , which (likewise) is below 
the Galactic average of 4.4-5.2 x 10 21 cm -2 mag -1 , but still within the tail of the distribution 
for Galactic sight lines (Welty et ah 2012). The corresponding ratio for the molecular gas, 
N(R 2 )/E(B — V) ~ 1.2 x 10 21 cm -2 mag -1 , is somewhat higher than the average Galactic 
ratio, but well within the Galactic distribution. The line of sight to NGC 5195 therefore 
seems to be probing interstellar material in M51 with a fairly high molecular content, but 
with a somewhat lower overall gas-to-dust ratio compared to that most commonly seen for 
sight lines in the Milky Way. 

The H i column density that we derive for the M51 material toward NGC 5195 (i.e., 
log iV(H i) rj 20.75) should be secure since we obtain the same result from both emission and 
absorption features. Normally, when one compares a column density deduced from 21 cm 
emission to that obtained from an absorption-line tracer, there is some ambiguity related 
to the exact location of the background source; it is not usually known how much of the 
emitting gas is being probed by the source used for the absorption-line measurements. In this 
case, however, we are confident that the nucleus of NGC 5195 probes the entire line of sight 
through the spiral arm of M51. Furthermore, there should not be any discrepancies between 
the H I column densities derived from the optical absorption and 21 cm emission features 
due to structure transverse to the line of sight since the beam size for THINGS (~6"; Walter 
et al. 2008) is comparable to the size of the slit used for the echellc observations (~3 , . , 2). 

It is therefore somewhat surprising that the H I 21 cm emission profile does not show 
better agreement in velocity space with the interstellar absorption features observed in the 
direction of NGC 5195 (Figure 3). The intensity-weighted mean H I velocity is ulsr — 
434 km s -1 , while the atomic and molecular absorption lines show a mean column-density 
weighted velocity of ulsr = 411 km s -1 . Figure 3 clearly demonstrates that the absorption 
peak is shifted to lower velocities compared to the peak in H I emission. Moreover, the 
extent of the H I emission seems to coincide better with the weaker absorption components 
seen to the red of the main components in the Na I, K I, Ca II, and CH + profiles. Since 
these latter species trace both atomic and molecular gas, it seems likely that the material 
“missing” from the H I profile at the velocity of the main interstellar absorption components 
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is mostly molecular. 

There are other indications that the line of sight to NGC 5195 is probing relatively dense 
gas in M51. Our lower limit on the Na I column density of log N(Na i) > 14.0, combined 
with our determination of log N(Ca n) = 12.96, yields N(Na l)/lV(Ca il) > 11. High 
N(Na l)/iV(Ca n) ratios such as this are usually associated with cold, dense quiescent clouds, 
where Ca is heavily depleted within interstellar grains. While the N(Na l)/W(Ca il) ratio can 
also be affected by variations in the ionization conditions (which can alter the abundances 
of both Na I and Ca il), the dominant effect in this case is probably the depletion of Ca. At 
log N (Na i) > 14.0, Ca II is likely to be the dominant ionization state of Ca in this material 
(e.g., Welty et al. 1996). We can therefore derive the Ca depletion factor directly from the 
Ca II column density by combining this value with our estimate for lV(H tot ). Assuming that 
the metallicity of the M51 material is roughly solar, and adopting the solar Ca abundance 
from Lodders (2003), we find [Ca/H] = log N (Ca il) — log W(H tot ) — log (Ca/H)© = —2.66± 
0.15, which is comparable to some of the largest depletions seen for Ti in the local Galactic 
ISM (e.g., Jenkins 2009; Welty & Crowther 2010). 

We have already remarked on the general weakness of the DIBs toward NGC 5195, which 
may be an indication that the gas being probed has a high molecular content. We further 
note that the ratio of the equivalent widths of the A5797.1 and A5780.5 DIBs is quite high; 
we find W\ (5797.1) /W\ (5780.5) = 1.55±0.26, which, to our knowledge, is among the highest 
values ever measured for this ratio (e.g., Vos et al. 2011; Welty et al. 2014). Such a high 
value for H / a(5797.1)/B / a(5780.5) may suggest either that the material is immersed in a very 
weak ambient radiation field, or else that the environment is significantly shielded from UV 
radiation (Vos et al. 2011). A low-UV environment for the M51 material could potentially 
explain why the total hydrogen column density predicted directly from N(K i) seems to be 
somewhat high (compared to our adopted value). A factor-of-two enhancement in N (K i), 
due to a weak radiation field, would be sufficient to explain the difference in the predictions 
for V(Htot)- As high W 7 a(5797.1)/W / a(5780.5) ratios also appear to be associated with high 
molecular fractions (Weselak et al. 2004), the high value of W'a(5797.1)/B / a(5780.5) that we 
find for the M51 material toward NGC 5195 is consistent (qualitatively) with our estimate 
of /(H 2 ) « 0.65. 

One complication to this picture is that the W(CH + )/W(CH) ratio toward NGC 5195 
seems also to be fairly high; our upper limit on the CH column density, combined with our 
measurement of log W(CH + ) = 13.91, gives W(CH + )/W(CH) > 2.3, while, for relatively 
dense and/or significantly shielded gas with a high molecular fraction, this ratio would be 
expected to be < 1.0 (e.g., Pan et al. 2005; Welty et al. 2006; Ritchey et al. 2006; Sheffer 
et al. 2008). A similar result was found for the ISM of M82 along the line of sight to 
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SN 2014J (Welty et al. 2014; Ritchey et al. 2015). The material in that direction was 
shown to have an unusual combination of a high H / a(5797.1)/H / a(5780.5) ratio and a high 
!V(CH + )/lV(CH) ratio, compared to sight lines in the Milky Way, in the LMC and SMC, and 
in other supernova host galaxies. Compared to the M82 gas toward SN 2014J, the material 
in M51 toward NGC 5195 may have a somewhat lower !V(CH + )/lV(CH) ratio and has a 
higher H7\(5797.1)/H / a(5780.5) ratio, but the M51 results still seem to be offset from the 
results for other sight lines (see, e.g., Figure 5 in Welty et al. 2014). 

The high CH + abundances in the M51 gas toward NGC 5195 and the M82 gas toward 
SN 2014J suggest that non-thermal chemical processes are highly active in these galaxies 
(e.g., Zsargo & Federman 2003). The formation of CH + occurs through the endothermic 
reaction C + + H 2 —» CH + + H, which cannot proceed at the low temperatures of diffuse 
molecular clouds (due to its activation energy of A E/k-Q = 4640 K). Thus, the large abun¬ 
dances of CH + observed in the Galactic ISM have traditionally been difficult to explain. 
Joulain et al. (1998) proposed that intermittent bursts of turbulent dissipation would be 
sufficient to heat the gas locally so as to account for the observed CH + abundances. Their 
model of turbulent dissipation regions (TDRs) has been further developed by Godard et 
al. (2009, 2014), who ford that the CH + abundance is directly proportional to the average 
turbulent dissipation rate. Tubulence is enhanced in the ISM of interacting galaxies (e.g., 
Irwin 1994; Elmcgreen et al. 1995), and both M51 and M82 are known to be interacting 
with their neighbors (M51 with NGC 5195 and M82 with M81). A possible explanation, 
therefore, of the high CH + abundances in M51 and M82 is that these abundances are linked 
to enhancements in interstellar turbulence resulting from galaxy-galaxy interactions. Such 
enhancements could be directly related to the gravitational encounters occurring between 
the interacting galaxies or could be due to subsequent increases in the star formation and 
supernova rates triggered by those interactions. 


5. CONCLUDING REMARKS 

We have discussed moderately-high resolution echellc observations of the nucleus of 
NGC 5195, the line of sight to which probes interstellar gas associated with the northeastern 
spiral arm of M51. Absorption from Na I, K I, Ca II, and CH + is detected at a velocity close 
to that exhibited by H I 21 cm emission from M51 at the position of NGC 5195. However, 
the H I emission profile (provided by THINGS) peaks at a slightly higher velocity, compared 
to the atomic and molecular absorption profiles, and coincides better with weaker absorption 
components seen to the red of the main components. Eight DIBs are detected at velocities 
similar to those of the atomic and molecular species seen in absorption. The H I column 
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density implied by the equivalent width of the A5780.5 DIB is consistent with that obtained 
from the THINGS integrated H I map of M51 at the location of NGC 5195. Intriguingly, 
the H 2 column density predicted from the observed column density of K I is comparable to 
IV(H i), suggesting a high molecular fraction (~0.65) for the M51 gas toward NGC 5195. 
The DIBs toward NGC 5195 are, on average, about 40% weaker than would be expected 
based on the K I column density, which is a further indication that the gas in this direction 
has a high molecular content. After correcting for the overall weakness of the DIBs, we 
estimate the reddening toward NGC 5195 to be E(B — V) 0.46 mag. The M51 material is 
characterized also by a high lV(Na i)/AT(Ca n) ratio (> 11), indicative of a high degree of Ca 
depletion, and a high W / a(5797.1)/H / a(5780.5) ratio (~1.6), suggestive of either a very weak 
ambient radiation held or a significantly shielded environment. A high N (CH + ) / N (CH) 
ratio (> 2.3) for the gas in this direction may be the result of enhanced turbulence due to 
interactions between M51 and its companion. 

While there have been other absorption-line studies of moderately dense gas in nearby 
galaxies, typically using a supernova as the continuum source (e.g., Cox & Patat 2008, 2014), 
an important aspect of this investigation is that the background source is not a supernova, 
and hence will not fade with time. This opens up the possibility of further high-resolution 
observations, which could explore in more detail the physical and chemical properties of the 
interstellar gas in M51 toward NGC 5195. For example, ground-based observations with 
higher S/N and higher resolution than those presented here could target absorption from 
Li I, and could potentially yield the first determination of the Li isotope ratio in a galaxy 
beyond the Local Group. A more comprehensive examination of the DIBs in this direction 
could yield important information concerning the behavior of the DIBs in an environment 
different from that typically encountered in the local Galactic ISM. Finally, observations of 
NGC 5195 with the Cosmic Origins Spectrograph onboard the Hubble Space Telescope would 
provide access to important UV absorption-line tracers, which could be used to study the 
abundances and physical conditions in the M51 material in much greater detail than can be 
done from the ground. Such information would be particularly useful in more confidently 
identifying the mechanism responsible for the high CH + abundance, and could ultimately 
provide valuable constraints for models of CH + formation in the ISM of our Galaxy and 
beyond. 


We thank Dan Welty for offering many useful comments and suggestions on this work. 
Support for this research was provided by the Kenilworth Fund of the New York Community 
Trust. The results presented here are based on observations obtained with the Apache Point 
Observatory 3.5-rn telescope, which is owned and operated by the Astrophysical Research 
Consortium. 
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Table 1. Absorption and Emission Features in NGC 5195 


Species 

A 

Aobs 

't ; LSR FWHM a 


(A) 

(A) 

(km s” 1 ) (km s -1 ) (km s -1 ) 


Stellar Absorption Lines 

Mg I 

5183.60 

5194.07 

+606 

252 

99 

Ca ii 

8498.02 

8515.80 

+627 

279 

118 


8542.09 

8559.58 

+614 

309 

128 


8662.14 

8679.95 

+616 

249 

103 

Emission Lines 

Ha 

6562.82 

6574.94 

+563 

157 

66 

[0 i] 

6300.30 

6312.45 

+578 

125 

53 

[N ii] 

6548.05 

6560.65 

+577 

171 

72 


6583.45 

6595.78 

+561 

194 

82 

[S ii] 

6716.44 

6729.16 

+568 

169 

72 


6730.82 

6743.66 

+572 

170 

72 
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Table 2. Absorption and Emission Features in the Nucleus of M51 


Species 

A 

-^obs 

■clsr FWHM a 


(A) 

(A) 

(km s -1 ) (km s _1 ) (km s -1 ) 


Stellar Absorption Lines 

Mg I 

5183.60 

5191.78 

+473 

206 

78 

Ca ii 

8498.02 

8511.71 

+483 

233 

98 


8542.09 

8555.89 

+484 

268 

110 


8662.14 

8675.70 

+469 

190 

77 

Emission Lines 

Ha 

6562.82 

6573.62 

+502 

182 

77 

H /3 

4861.33 

4869.65 

+513 

108 

46 

[0 i] 

6300.30 

6310.64 

+492 

138 

59 

[0 III] 

4958.91 

4967.38 

+512 

166 

71 

[0 m] 

5006.84 

5015.23 

+502 

195 

83 

[N II] 

6548.05 

6558.96 

+500 

184 

78 


6583.45 

6594.41 

+499 

216 

92 

[S ii] 

6716.44 

6727.57 

+497 

176 

75 


6730.82 

6742.10 

+503 

172 

73 









Table 3. 


Interstellar Atomic and Molecular Absorption Lines toward NGC 5195 


Species 

A a 

(A) 

log /A a 

(^LSR.) b 
(km s _1 ) 

w x 

(rnA) 

log N a 

log -ATadopt 

Ca ii 

3933.661 

3.392 

+405 

371 ± 39 

12.73 + 0.05 

12.96 + 0.08 


3968.467 

3.092 

+410 

242 ± 43 

12.84 + 0.07 


Na i 

5889.951 

3.577 

+415 

1556 ± 8 

13.27 + 0.01 

>14-0 


5895.924 

3.276 

+414 

1403 ± 8 

13.49 + 0.01 


K i 

7664.899 

3.710 

+413 

325 + 6 

12.03 + 0.01 

12.11 + 0.01 


7698.965 

3.409 

+413 

191 + 5 

12.07 + 0.01 


CH+ 

4232.548 

1.363 

+410 

67 + 16 

13.91 + 0.10 

13.91 + 0.10 

CH 

4300.313 

1.338 


<30 

<13-6 

<13-6 


a Wavelengths and /-values are from Morton (2003) for atomic lines and Gredel et 
al. (1993) for molecular lines. 

b Column-density weighted mean velocity: (ulsr) = x v)/( ^ N a (v)). 

c Adopted column density, including any correction for line saturation. 







Table 4. Diffuse Interstellar Bands toward NGC 5195 


A 

(A) 

\ a 
^obs 

(A) 

VLSR 
(km s -1 ) 

FWHM a 

(A) 

fb+meas 

(mA) 

II+pretD 

(mA) 

E(B-V) pred d 

(mag) 

E(B — V) COII e 
(mag) 

5780.5 

5788.6 

+418 

1.49 

71 ± 10 

288 

0.13 

0.36 

5797.1 

5805.1 

+414 

1.07 

110 + 10 

121 

0.56 

0.54 

5849.8 

5857.8 

+408 

1.02 

42 + 6 




6196.0 

6204.4 

+408 

0.28 

10 + 2 

33 

0.16 

0.35 

6203.6 

6212.2 

+416 

3.87 

126 ± 16 

110 

0.60 

0.55 

6376.1 

6384.7 

+406 

1.11 

37 + 6 




6379.3 

6388.0 

+407 

0.56 

26 + 5 

55 

0.29 

0.49 

6613.6 

6622.7 

+414 

1.16 

75 + 6 

133 

0.38 

0.46 


a Observed central wavelengths and band widths of the DIBs were determined according 
to the definitions given in Hobbs et al. (2008). 

b Measured Wa(DIB). 

c Predicted Ha(DIB), assuming log IV(K i) ~ 12.11 and using Galactic relationships from 
Welty et al. (2014). 

d Predicted E(B — V ) from measured Ha (DIB) and using Galactic relationships given in 
Welty et al. (2014). 

e Predicted E(B — V), adjusted for residuals, using Galactic relationships given in Welty 
et al. (2014). 
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Fig. 1.— Final co-added spectrum of NGC 5195 in the vicinity of Na I D. Narrow Galactic 
Na I absorption components are seen near the rest wavelengths of the D lines. The broader, 
redshifted absorption features are presumably due to interstellar material in the spiral arm 
of M51. These features are superimposed onto very broad stellar Na I D lines, which are 
further redshifted and presumably arise from red giant stars in the nucleus of NGC 5195. 
The dashed line shows the continuum fit adopted for the purpose of analyzing the interstellar 
features. Strong Na I night sky emission lines are identified with the symbol ©. 





Fig. 2.— Final co-added spectrum of NGC 5195 in the vicinity of the (redshifted) A5780.5 
and A5797.1 DIBs (which can be seen at wavelengths of 5788.6 and 5805.1 A, respectively). 
The dashed line shows the adopted continuum obtained through a low-order polynomial 
fit to the observed spectrum (excluding the regions showing DIB absorption). Also shown 
(displaced by +0.1 in relative intensity units) is a smoothed spectrum of the prototypical red 
giant star Arcturus (Hinkle et al. 2000), which provides a template for the underlying stellar 
absorption features in the NGC 5195 spectrum. The Arcturus spectrum has been shifted by 
+616 km s -1 in accordance with the mean velocity exhibited by the stellar absorption lines 
in NGC 5195 (see Section 3.1). 
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Fig. 3.— Continuum-normalized interstellar absorption profiles of Na I, K I, Ca II, and 
CH + in the direction of NGC 5195, along with the H I 21 cm emission profile of M51 at 
the position of NGC 5195 from THINGS. Spectra covering the CH region is also shown. 
Dotted features in the bottom panel denote overlap in the Na I D lines. The vertical dotted 
line marks the mean column-density weighted velocity of the interstellar species seen in 
absorption (ulsr = 411 km s -1 ). The mean H I velocity in this direction from the THINGS 
velocity field (moment 1 map) for M51 is ulsr = 434 km s -1 . 
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Fig. 4.— Continuum-normalized absorption profiles of the eight DIBs detected in the di¬ 
rection of NGC 5195. The vertical dotted line marks the mean central velocity of the DIBs 
(vlsr = 411 km s -1 ). Solid points indicate the limits of integration used to determine the 
equivalent widths. Note that the asymmetric profile of the A6203.6 DIB, which includes an 
extended red wing, closely matches the profile shape of this DIB in Galactic sight lines (e.g., 
Friedman et al. 2011). 
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Fig. 5.— Continuum-normalized emission profiles for prominent emission lines seen in the 
spectrum of NGC 5195. The continuum region surrounding the [N n] and Hct emission 
lines is affected by broad stellar Ha absorption, which had to be fitted and removed from 
the spectrum before proceeding with the normalization. The vertical dotted line marks the 
mean centroid velocity of the emission lines listed in Table 1 (ulsr — 570 km s' 1 ). 
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Fig. 6. — Continuum-normalized absorption profiles for the Ca II infrared triplet arising 
(presumably) from red giant stars in the nucleus of NGC 5195. The vertical dotted line 
marks the mean centroid velocity of the stellar absorption lines listed in Table 1 (ulsr = 
616 km s -1 ). 








